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The systematic treatment of heavy quark mass effects in DIS in current CTEQ global 
analysis is summarized. Applications of this treatment to the comparison between 
theory and experimental data on DIS charm production are described. The possibility 
of intrinsic charm in the nucleon is studied. The issue of determining the charm mass 
in global analysis is discussed. 

1 Introduction 



Contemporary global QCD analyses of high precision Deep Inelastic Scattering (DIS) data, 
along with other hard processes, require a consistent treatment of heavy quark mass effects 
in the perturbative QCD (PQCD) framework. This review pQ summarizes key features of 
the formalism implemented in the current CTEQ global analysis project [3] and results on 
its application to heavy flavor physics in global analysis j2j. Sec.[3]presents the results of the 
^vq ' new global fits compared to heavy flavor production data in DIS [3j. Sec. H] addresses issues 

(f) . related to possible intrinsic charm in the nucleon [BJ. Sec. H] discusses the topical question: 

' can the charm mass be reliably determined in global QCD analysis? 

[~^. . Due to space limitation, it is impossible to include in this short written report the 

figures that illustrate the results discussed in the corresponding talk, as summarized above. 
However, since the slides for the talk have been made available at the official conference URL 
PQ , we shall make use of these, and refer the reader to the actual figures by the slide numbers 
where they appear in the posted talk pQ. The same space limitation restricts citations to 
only the papers and talks on which this report is directly based. 
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2 General PQCD framework including heavy quark masses 

The key features of the general-mass PQCD framework of [3] is illustrated in slide 3 of [T]. 
Factorization Formula and (scheme-dependent) summation over parton flavors: 

Collins has established that the PQCD factorization theorem for the structure functions 
takes the general form F\(x, Q 2 ) = Yl a f a ® even when the heavy quark mass effects 
are kept. Here, the summation is over the active parton flavor label a, f a (x,/j,) are the 
parton distributions at the factorization scale /i, and (x, Mi/ ix] are the infrared safe 
Wilson coefficients (or hard-scattering amplitudes) that can be calculated order-by-order 
in perturbation theory. The summation over "parton flavor" label a in the factorization 
formula is determined by the factorization scheme chosen to define the parton distributions 
/ a (x,/i). In general, we use the variable flavor number scheme. 

The summation over (physical) final-state flavors: For total inclusive structure func- 
tions, the factorization formula contains an implicit summation over all possible quark flavors 



DIS 2007 



in the final state: u> a = J^b^ai where "&" denotes final state flavors, and uj\ is the pertur- 
batively calculable hard cross section for an incoming parton "a" to produce a final state 
containing flavor "6 " . It is important to emphasize that "6" labels quark flavors that can 
be produced physically in the final state; it is not a parton label in the sense of initial- 
state parton flavors described in the previous subsection. In a proper implementation of 
the general-mass (GM) formalism, the distinction between the initial-state and final-state 
summations must be unambiguously and correctly observed. 

Kinematic constraints and rescaling: Kinematic constraints from the phase space treat- 
ment have a significant impact on the numerical results of the calculation. In DIS, with heavy 
flavor produced in the final state, the most natural way to ensure the correct kinematics for 
both NC and CC processes is to use the rescaling variable \ = x(l + (S/ Mf/Q) 2 ) in place 
of the usual Bjorken x in the convolution integral of the factorization formula. Here Sj Mf 
is the sum of all heavy flavor masses in the final state. This is the ACOT% prescription used 
in most recent literature. 

Hard Scattering Amplitudes and the SACOT Scheme: The hard scattering ampli- 
tude (x, Q/n, Mi/[i) is by definition infrared safe, meaning it is free from logarithmic 
"mass-singularities" in the limit M,/Q — > 0. Within the PQCD formalism, there is some 
freedom to choose how the finite mass effects are treated. The choice that makes the cal- 
culation simplest while retaining full accuracy (the SACOT scheme) can be stated as: keep 
the heavy quark mass dependence in the Wilson coefficients for partonic subprocesses with 
only light initial state partons (g, u, d, s); but use the zero-mass Wilson coefficients for sub- 
processes that have an initial state heavy quark (c, b). For the 4-flavor scheme to order a s 
(NLO), we do the following: (a) keep the full M c dependence of the gluon fusion subprocess; 
(b) for NC scattering (j/Z exchanges), set all quark masses to zero in the quark-initiated 
subprocesses; and (c) for CC scattering (W^ exchange), set the initial-state quark masses 
to zero, but keep the final-state quark masses on shell. 

Choice of Factorization Scale: The total inclusive structure function Ff ot is infrared 
safe. Consider the simple case of just one effective heavy flavor charm, F\ ot = F^ 9 * + F? 
for any given flavor-number scheme. Since the right-hand side of this equation is dominated 
by the light-flavor term F^ lght ; anc j the natural choice of scale for this term is yU = Q, it is 
reasonable to use this choice for both terms to ensure infrared safety. On the other hand, in 
the case of experimentally measured semi-inclusive DIS structure functions for producing a 
charm particle in the final state, F£ is theoretically infrared unsafe beyond NLO. One may 
nonetheless perform comparison of NLO theory with experiment with the understanding 
that the results are intrinsically less reliable, and they can be sensitive to the choice of 
parameters. The most natural choice of factorization scale in this case is [i = \JQ 2 + M 2 . 

3 Results and Comparison with heavy flavor production data 

Slides 4 and 5 of [I] show the size of heavy quark mass effects on the calculation of i^Oc, Q) 
and Fl(x, Q). The color coded areas (with complementary contours) indicate the fractional 
differences between GM and zero-mass (ZM) calculations. Understandably, the largest dif- 
ferences occur at low Q and low cc; and the significance is much more for Fl(x,Q) than 
for F2(x,Q), since the former vanishes at LO for the ZM case. As indicated in slide 6, 
the GM calculation is stable and robust. It has been used as the basis for a new round of 
global analysis of PDFs, using the full set of HERA Run I neutral current (NC) and charged 
current (CC) total cross section and heavy flavor production data, along with the usual DY 
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and jet data cf. [3J- Here we shall only present the comparison of the new fits to the heavy 
flavor production data measured at HERA. 

Slide 10 shows the comparison of the ZEUS f996-97 and 1998-2000 charm production 
data to the theory values obtained with the new PDF sets CTEQ6.5M (same shape for 
strange and non-strange seas, [3]), CTEQ6.5S0 (independent shapes for strange and non- 
strange seas, [5]) as well as for the older CTEQ6HQ. Plotted are ratios of F-f (x, Q) to that of 
a best fit to the respective data set. The fits to data are all reasonable. The new PDFs give 
slightly better fits than the previous one. Slide 11 shows the comparison of the HI charm 
and bottom production data to the theory values from the same PDF sets. The F£(x,Q) 
data points have more scatter around the (smooth) theory values. The overall % 2 of these 
fits is however acceptable. 

It is worth noting that correlated systematical errors are always taken into account in 
our global analysis. The data points shown on these plots have been shifted by the fitted 
systematic errors; hence the differences between the data points and the theory values as 
they appear on these plots give a faithful indication of the quality of the fits. 

4 Is there intrinsic charm in the nucleon? 

Many nonperturbative models of nucleon structure suggest the existence of intrinsic charm 
(IC) — a non-vanishing component of nucleon parton structure at the scale of M c . On the 
other hand, practically all global analysis of the parton structure of the nucleon so far 
ignore this possibility and make the simplifying assumption that all heavy quark partons are 
radiatively generated: they only arise from perturbativc QCD evolution, starting from zero 
at fi ~ M c . Where does the truth lie? The resolution of this dichotomy is of inherent physics 
interest because it concerns the fundamental structure of matter, as well as of practical 
interest because the cross sections for many beyond-the-standard-model (BSM) processes at 
hadron colliders depend on the charm parton content of the nucleon. We have addressed this 
problem phenomcnologically by a careful global analysis based on the GM PQCD formalism 
that, for the first time, allows for an independent charm sector [B]. 

As indicated in slide 14 of pQ, the following specific scenarios for the charm sea, c(x, /j, = 
M c ), are explored within our GM global analysis framework: (i) the conventional radiatively 
generated charm; (ii) non-vanishing IC c(x, M c ) that is sea-like (i.e. shaped as the light sea 
quarks); and (iii) IC of the kind suggested by light-cone wave function models of the nucleon 
(peaked at moderately large x). Within scenario (iii), we further distinguish two models: 
the one studied by Brodsky et al. (the BHPS model), and a meson cloud model. 

Slide 15 summarizes the main results. The figure shows the goodncss-of-fit for the global 
analysis, Xgiobal' as a function of the magnitude of the IC component, measured by the 
momentum fraction carried (x) c+s , under the various scenarios. In the range < (a;}c+c < 
0.01 (outlined by the horizontal oval), Xgiobal ^ s largely insensitive to (x) c +ci indicating that 
there is no strong evidence for or against IC of a magnitude in this range. However, outside 
this range, for (x) c+s > 0.01 (outlined by the vertical oval), we see a precipitous rise of X 2 g i bai 
as (x) c+E increases. Thus our global analysis sets a useful upper bound on the amount of 
intrinsic charm that is consistent with existing data. Using a 90% confidence level (C.L.) 
criterion, this bound is (x) c +c < 0.02. 

Although models of IC generally do not predict (x) c + E , typical guesstimates place it 
around 0.01. This is consistent with the bound we determined from the above global analysis. 
The presence of IC of such a magnitude can have an impact on certain BSM processes, such 
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as charged Higgs production in hadron collider phenomenology. Cf . slide 7 of [5] , and [H [5] . 
Slide 16 shows the charm distribution c(x, fx) at three energy scales /i = 1.3, 3.16, 85 GeV in 
the BHPS scenario, for various magnitudes of the initial distribution. We see the radiatively 
generated component (peaked at small x) catching up with the IC component (peaked at 
moderate x) as \i increases. However, the latter clearly still dominates in the x region ^0.1 
even at the W/Z mass scale. 

5 Can the charm mass be determined in global analysis? 

In principle, heavy quark masses Af 2 ;(/i) at some renormalization scale /i are basic parameters 
of QCD, similar to the coupling a s (p,). Thus, just like for a s , there has been recent interest 
in determining Mj, e.g. the charm mass M c , from global QCD analysis. In particular, is 
it possible to perform a conventional global QCD analysis using M c as one of the fitting 
parameters, and thereby determine the charm mass to be the one that gives the best fit? If 
so, one may further ask, is this mass the MS-mass or the pole-mass? 

Slides 18-20 show results of a study, following the above procedure literally: one finds 
that the global analysis favors a relatively small values of M c ^1.3 GeV, and the goodncss- 
of-fit Xgiobal increases with M c . But, a closer examination of the problem immediately raises 
the question: what is the physical meaning of this favored value of M c ? The problem is, 
a chosen value of M c affects the global QCD analysis in two distinct ways: (i) through 
the mass-dependent Wilson coefficients in the theoretical calculation (the pole-mass); and 
(ii) through the initial condition c(x, /i = M c ) = — the implicit assumption of radiatively 
generated charm that is used in all existing global analyses. It turns out, the global fit is 
influenced much more by the latter than by the former. Since radiatively generated charm 
is only an assumption, not an integral part of the QCD theory, the value of M c favored by 
global analysis is not directly related to the basic QCD charm mass parameter — it is neither 
the M^-mass nor the pole-mass! In order to answer the original question "can the charm 
mass be determined in global analysis?" , one needs to clearly differentiate between the two 
sources of dependence on M c mentioned above. This is currently under study. 
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